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^ Abstract 

Coherent transport of atoms trapped in an optical lattice can be controlled by microwave-induced 
spin flips that correlate with site-to-site hopping. We study the controllability of homogeneous 
one-dimensional systems of noninteracting atoms in the absence of site addressability. Given these 
restrictions, we construct a deterministic protocol to map an initially localized Wannier state to a 

Oh 

wave packet that that is coherently distributed over n sites. This is extended to analytic solutions 
for arbitrary unitary maps given homogenous systems and in the presence of time-dependent uni- 
form forces. Such control is important for applications in quantum information processing such as 
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quantum computing and quantum simulations of condensed matter phenomena. 
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I. INTRODUCTION 



Neutral atoms trapped in optical lattices have emerged as a rich platform for exploring 
a wide variety of phenomena and devices based on coherent quantum dynamics. Examples 
include quantum computers [1, 2, 3, 4], quantum simulators of condensed matter [5, 6], topo- 
logical quantum field theory [7, 8], and quantum chaotic dynamics [9, 10, 11]. An essential 
ingredient in these systems is the coherent control of atomic transport in the lattice. Such 
transport is driven by time-dependent variations in the lattice potential and the application 
of external fields. In its most basic form, the atoms' ballistic tunneling between sites in 
a sinusoidal potential can be controlled through time-dependent modulations of the lattice 
depth and phase. The latter can be used to impart a time-dependent acceleration to the 
lattice, thereby simulating the effects of an applied electric field for electrons in a crystal that 
give rise to the fundamental paradigms of coherent transport in solid-state physics. Bloch 
oscillations [12], Wannier-Stark ladders [13], Landau-Zener tunneling [14], and dynamical 
localization [15] have all been demonstrated in optical lattices and explored as mechanisms 
for coherent control. 

More complex lattice geometries introduce additional features. For example, in a lattice 
of double- wells, one can drive transport between sites in a pairwise manner, assuming a 
sufficient barrier to ignore tunneling between different double- wells [16, 17, 18, 19]. In 
this case, the control problem is substantially simplified, as the relevant Hilbert space in 
a given time interval is restricted to a small discrete set of energy levels, as opposed to 
the infinite chain of levels in a sinusoidal lattice. Control across the entire lattice can be 
implemented by modifying the geometry so that the wells are alternatively coupled to all 
nearest neighbors (left or right in ID). Designer double- well lattices have been explored for 
quantum information processing tasks such as quantum computing [20] and simulations of 
condensed-matter phenomena [21]. 

Still richer control is possible for spinor lattices where the optical potential depends on 
the atom's internal spin state [22]. The lattice's morphology can now be modified through 
variation of the laser polarization as well as intensity, lattice phase, etc. The earliest pro- 
posals for quantum logic in optical lattices via controlled collisions involved transport of the 
atoms via time-dependent rotation of the direction of a laser beam's polarization [1, 2, 23]. 
An alternative and perhaps more robust route to coherent control of atomic transport is to 
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use external fields to drive spin- changing transitions that are correlated with atomic motion, 
similar to the scheme proposed by Foot et al [24] . Such protocols can make use of the tools 
for robust control of spins [25, 26, 27], as developed in NMR, to the control of atomic motion 
in the lattice. 

In this article we explore methods for coherent control of atomic transport with 
microwave-induced spin rotations between hyperfine levels and polarization-gradient lat- 
tices. Our main focus is on controllability - how the Hamiltonian that governs the dynamics 
restricts the possible unitary maps that one can implement, and how to design specific wave- 
forms to carry out a given task. We will consider here the simplest problem of noninteracting 
atoms in one dimension. While extensions to the interacting case are nontrivial, the current 
work is an important stepping-stone in that direction. 

The remainder of the paper is organized as follows. In Sec. II we establish the formalism 
necessary to describe spinor lattices and their interaction with external fields. We apply 
this to study the conditions for wave function control (the preparation of a desired spinor 
wave function starting from a known localized Wannier state) and prescribe a constructive 
algorithm for carrying out this task in Sec. III. We then generalize this in Sec. IV to the 
case of more general unitary maps for unknown initial states. Finally, we summarize and 
give an outlook towards future research in this area in Sec. V. 

II. MICROWAVE-DRIVEN SPINOR LATTICES 

Spinor lattices arise from the tensor nature of atom-photon interaction. In a monochro- 
matic laser field Re (E(x)e z ^^ _a;L ^) , the light-shift potential takes the form, 

V LS (x) = -\a l3 E:{s)E 3 {s), (1) 

where is the atomic dynamic polarizability at frequency uol, for atoms in a particular 
ground-state manifold. We consider a ID geometry consisting of counterpropagating laser 
beams with linear polarizations, forming a relative angle 9 (the "lin-^-lin" geometry). Re- 
stricting our attention to alkali atoms, when the laser field detuning is large compared to 
the excited-state hyperfine splitting but not large compared to the fine-structure splitting, 
ellipticity in the laser field leads to a fictitious magnetic field that varies periodically in 
space, and results in a spin-dependent light shift [28]. Taking the quantization axis along 
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one of the laser beams' wave vectors, the spinor light-shift potential can then be written 

V LS = ^2[E f - V F>mF cos {2k L z + S Fmp )} \F, m F ) (F, m F \ , (2) 

where E F is the degenerate energy of a hyperfine manifold, Vp ?mF = 
Vo^/cos 2 9 + ((/f^f/2) 2 sin 2 0, and S^ mF = tan -1 [^f^f tan 0/2). Here, Vo is the lat- 
tice depth for lin-||-lin and gp is the Lande g-factor of the hyperfine manifold under 
consideration. 

The addition of a static bias magnetic field breaks the degeneracy between Zeeman sub- 
levels within a manifold and allows us to spectrally isolate different microwave transitions 
between the manifolds. For concreteness, we consider 133 Cs, and choose a two-state sub- 
space \F — 4, to = 3) = ||) , \F — 3, to = 3) = ||) to define a pseudo-spin- 1/2 particle. Re- 
stricting to this subspace, and adding a near-resonant microwave field with magnetic field 
cos (uo^wt — (j)) that couples these spin states, the Hamiltonian in the rotating frame 
takes the form H = Hi att + H^ W) where 

2 

Hiatt = £ V cos (2k L z + 5 ) I TXT I - V cos (2k L z - 5 ) | |> (| I (3a) 

H^ w = ^cr z (cos (f)^ w a x + sin (f)^ w a y ) . (3b) 

Here Vo = ~ V^, S = « —£3,3, = (T| A ID B^ w /fi is the microwave resonant 
Rabi frequency, is the microwave detuning from a hyperfine resonance defined by the 
untrapped atoms, c/)^ is the phase of the microwave oscillator, and the Pauli-a operators 
are defined relative to the pseudo-spin. 

Neglecting the kinetic energy and diagonalizing the Hamiltonian leads to adiabatic or 
microwave- dressed potentials , 

V±(z) = -V cosS cos(2kz) ± ^y/(2V sin S sm(2kz) - A^) 2 + n% w . (4) 

At = 0, in the lin_Llin (0 = tt/2) geometry, the adiabatic potentials yield a period A/4 
"subwavelength" lattice [29, 30, 31]. In the context of a Hubbard Hamiltonian describing 
interacting particles moving on a lattice [32] , this configuration gives us greater freedom to 
independently control the site-to-site tunneling rate J and the onsite interaction strength U 
[33]. By employing both optical and microwave fields, the lattice depth dominates the control 
of U while the the applied microwave dominates control of J. Moreover, the tunneling matrix 
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element is complex, set by the microwave phase, allowing for time-reversible tunneling and 
further control [34]. 

For 9 ^ mr/2, the adiabatic potentials take the form of a lattice of double- well potentials 
arising from the asymmetry for transport to the left vs. the right. The parameters char- 
acterizing the double well, including barrier height, tunneling matrix element, and energy 
asymmetry ("tilt"), can be controlled through variations of lattice intensity/polarization, 
microwave power, and detuning. The richness of this system should enable us to control 
wave function coherence for spinors over multiple sites. Our early work on this subject 
demonstrated spinor double- well coherence driven by Larmor precision in a quasistatic mag- 
netic field [22]. The current approach, based on applied microwave fields should be much 
more robust and controllable. 

While the dressed-lattice adiabatic potentials guide intuition about the transport, quan- 
titative predictions are more accurately made by considering the band structure of the 
Hamiltonian, Hi att) in Eq. (3). Associated with the spin s =| and s —[ lattices are Bloch 
states for band-n and quasimomentum-g, \^n\)-> and Wannier states for that band and lattice 
site-/, related by the usual Fourier transform over the first Brillouin zone, 



Here and throughout, lengths are measured in units of the lattice period L = X L /2 and wave 
numbers in units of the reciprocal lattice vector K = Att/Xl. For sufficiently deep lattices 
and atoms in the lowest lying bands, tunneling between sites is completely negligible over 
the timescales of interest. In that case, the lattice Hamiltonian is diagonal both in the Bloch 
and Wannier bases, with no energy variation over the q or I index. 

Transport dynamics are driven by the microwaves tuned to cause transitions between 
the ground bands associated with the spin-up and spin-down lattices. We assume that the 
detuning and Rabi frequency are sufficiently small so that single-band/lattice tight-binding 
(TB) model is a good approximation. Henceforth we drop the band index and set n — 0. In 
the Wannier basis, the total Hamiltonian in the TB approximation is 
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(5) 
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Htb = ^ 
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where, 



oi = |i,T) </,T|-|U> <UI, 



(7a) 
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i,T> <UI 



(7c) 



are the Pauli operators for two-level transitions that pairwise couple spin-down Wannier 



Note, we have chosen an arbitrary labeling of the Wannier state indices by convention so 
that a spin-down state and spin-up state to its right are both associated with the same 
lattice period label, I. Because the microwaves transfer negligible momentum to the atoms, 
translation of the atomic wavepacket is possible only when the probability amplitude of an 
atom overlaps between neighboring sites. The Rabi frequencies for transitions to the left or 
right are thus weighted by Franck-Condon factors, Qr = (^Jl^jf)^^ = (^j-il^i)^^- 
For the ground bands in the TB approximation, a large asymmetry in right-left transport 
and isolation of double wells arises from small asymmetry in right-left displacement of the 
lattice due to the Gaussian overlap of the wavepackets (see Fig. 1). 

The combination of spinor optical lattices and microwave- driven spin dynamics provides 
a wide variety of parameters that can be modulated in real time during an experiment to 
coherently control atomic transport. In the next section we study the formal controllability 
of this system and develop constructive protocols to implement desired unitary maps. 

III. WAVEPACKET CONTROL 

Given the time-dependent Hamiltonian at hand, a first question to address is "controlla- 
bility", i.e., which class of unitary transformations can be generated by the arbitrary design 
of the waveforms that parameterize that Hamiltonian. We consider first the problem of 
preparing an arbitrary wavepacket that is coherently distributed over multiple sites of the 
lattice, starting from an initially localized Wannier state. Unless specially designed to allow 
for individual site addressability [4, 35], control in a typical lattice is limited by translational 
invariance of the operations. Although a linear gradient breaks the translational symmetry, 
the controllability of the system is still limited, as we show below. We thus restrict our 
attention to strictly periodic lattices and allow for an additional constant force F on the 
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6 = 90" = 80" 




FIG. 1: Controlled transport in lin-0— lin. Spin-dependent lattices have a relative phase shift 
due to polarization gradient. In the lin±lin configuration, 9 — 90°, there is no asymmetry in 
right-left transport, and atoms can ballistically tunnel in the dressed potential. For a change of 
just 10° degrees away from lin±lin, in a lattice with an oscillation frequency of 20Er, the ratio of 
the effective tunelling rates is Ql/^r ~ 3500. The result is lattice of double wells with pair- wise 
tunneling couplings. 

atoms, as in the case of lattices held vertically in gravity, or when the overall lattice is 
accelerated through time-dependent changes of the standing wave pattern. 

We consider Hamiltonians which are composed of a translationally invariant part, 
with period L, and an applied force, F, 

H(t) = H (t) + F(t)x (8) 

A particular example is H (t) = H TB (t), as given in Eq. (6), with time-dependent variations 
in microwave power and/or phase. The time evolution of such a system may be written in 
the form of the time-ordered exponential 

U(t) =Texpj-i^ (H Q (t') + F(t')x)dt'^, (9) 

Such a map has the property that if we translate the entire system by j times the period of 
H , then U e'^o F(«)dtjL Um If the initial state i^\ 

is a localized Wannier state, it satisfies 
WW) = <5 il0 , (10) 
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where 7} translates the system by jL. If state \(j>) maps to \(j)') under a unitary evolution of 
this form, then the evolved state satisfies this same condition, as follows from the identity 

(<j>'\T0) = (d>\^TjU\<f>) = e i X F W L (<l>\T j \<l>) = S jt0 . (11) 

Thus any wavepacket prepared by these controls must be orthogonal to itself after a trans- 
lation by an arbitrary number of periods. Furthermore, since the force F has dropped out, 
the linear gradient does not impact the range of states that my be reached. 

Are all states that satisfy this constraint reachable through some choice of the control 
waveforms that parameterize Htb i n Eq. (6)? To show that this is the case, we employ a 
protocol for constructing a desired state-to-state mapping as defined by Eberly and Law in 
the context of Jaynes-Cummings ladder [36]. First note that if we can map the state \(j>) to 
another state and the control Hamiltonian allows the unitary map to be time-reversible, 
then we can map \(j) f ) to \(j>). Thus, in order to show that we can get from a localized state 
to any state satisfying Eq. (11), we consider the time-reversed problem of mapping such a 
state to the initial Wannier state. 

To construct the desired map, we employ a series of SU(2) rotations on resolvable sub- 
spaces of the total Hilbert space. Such a collection of disjoint two-level systems can be 
addressed in a lin-0-lin spinor lattice in either an asymmetric configuration (0 ^ vr/2), or 
in a lin_Llin configuration in the presence a sufficiently strong uniform force so that isolated 
pairs of states are spectroscopically addressable (see Fig. 2). Note that in the latter case, 
because of the presence of an external force, if we translate the system by an integer mul- 
tiple of A/2, the state picks up an extra phase due to the linear gradient. As we will see 
below, our construction does not require these phases, and since the construction is capable 
of synthesizing all reachable states, the phases are redundant. We can ignore the phases if 
we choose the time over which the two level unitaries operate to be an integer multiple of 
2n/FL, and we will assume this to be the case for the rest of this section. 

We will restrict our attention to states with support strictly on a finite set of lattice sites 
and zero probability amplitude outside some range. Such spinors can be represented as 

m= E (qiIU)+q t I't»- (12) 

l—lmin 

The finite extent of the wave function, together with the condition expressed in Eq. (11), 
places a constraint on the two-level subspaces. If we translate the entire state by l max — l min 
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FIG. 2: Two different methods for isolating unitary maps on two-level subspaces. In (a) and (b), 
the choice of polarization angles in a lin-#-lin lattice isolate different sets of two-level systems with 
transport either to the right or to the left. In (c) and (d), two-level systems in a lattice with a 
linear gradient are spectrally addressed through their distinct microwave transition frequencies, 
which differ by ±FL/2Ti 



then, 



T, 



(13) 



Thus, in order to satisfy Eq. (11), the two outermost two-level subspaces must be orthogonal. 
Moreover, because of translational invariance, the unitary transformations that we apply are 
equivalent at each period of the lattice. In particular, by unitarity, if we apply a rotation 
operator that maps the subspace on the left end of the atomic distribution to pure spin up, 
the subspace at the right end of the distribution must be rotated to pure spin down. 

These observations are the core of our construction (see Fig. 3). A sequence of two-level 
rotations can be used to map a coherent superposition delocalized across the lattice to one 
localized at a single site in a single spin state. In the first step, a rotation is applied to 
map all population at the left-most two-level system (/ = l m i n ) to spin-up according to the 
microwave- driven SU (2) transformation, 
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1 



Ci + Cr I 



(14) 



where N = \q 



Q 



\ By the translation symmetry this simultaneously maps 



Initial State Step 1 




FIG. 3: An example of state preparation through sequential SU{2) rotations. The polarization 
of the two counterpropagating beams is chosen to isolate sets of two-level systems. In step 1, 
population in the leftmost two-level system is mapped to entirely spin up, forcing population in 
the right most two level system to be entirely spin down by Eq. (13). In step 2, the polarization 
is rotated so that a different set of two level systems are coupled. In step 3, the population in the 
leftmost level is then mapped to entirely spin down. In step 4 the polarization is set to its original 
configuration and steps 1-4 are repeated until the state is localized. For this particular case, we 
have chosen to end the sequence spin down. A different choice for the final pulse would have ended 
the state spin up. 

all population at the right-most two-level system (Z = l max ) to spin-down, resulting in the 
spinor state 

lmax 1 

1^') = C L„T \ l min T) + £ (4 i \ll)+4 l \ll))+<i na . i \La*l), (15) 

thereby shrinking the extent of the wave function by A/2. The lattice is then reconfigured 

10 



(either through polarization rotation, a change of microwave frequency, or a change in ac- 
celeration) so that the opposite neighbors are coupled (spin-down now coupled to spin-up 
neighbor on left). An appropriate 577(2) rotation is then applied to map all population on 
the left most edge to spin down (simultaneously moving all population on right-most edge 
to spin-up). Repeating, we form a sequence of rotations that take the outer edges of the 
distribution and map them inwards in steps of A/4 until all the population is localized at 
one Wannier state. Reversing the order of the sequence provides the desired protocol for 
constructing any single particle wave function in the given band and with support on a finite 
number of lattice sites, subject to the constraint Eq. (11). 



IV. IMPLEMENTING GENERAL UNITARY TRANSFORMATIONS 

In the previous section we studied the construction of a particular unitary transformation 
— mapping an initially localized Wannier state to a spinor wave function delocalized over a 
finite number of lattice sites, under the constraint of Eq. (9). In this section we consider the 
most general unitary map that we can implement under these constraints. We begin with 
the case of perfect translational invariance of the lattice. In addition, since the microwave 
photons possess negligible momentum, only Bloch states with the same quasimomentum are 
coupled. Because of these symmetries, any unitarity transformation that we can synthesize 
will be block diagonal in the Bloch basis, U = ® q U q) where each block is a U(2) matrix 
that connects spin-up and spin-down states with the same quasimomentum q. In the basis 

{|U),k,!>}, 

r a(q) f3*(q) 
J{q) -a*(q) 

where |^(g)| 2 + |/3(g)| 2 =l. If the two lattices are sufficiently deep, tunneling is suppressed 
so 7 is independent of q and can be factored out of the problem, leading to SU (2) rotations 
in each block. 

The decomposition of a translationally invariant unitary transformation into a direct sum 
of SU (2) matrices has important implications for the design of arbitrary maps. Generally the 
design of a time-dependent waveform that generates an arbitrary unitary map is substantially 
more complex than a protocol for state-to-state mapping on an initially known state [37]. 
Intuitively, this is because state-to-state maps only constrain one column of a unitary matrix, 



U q = e l ^ q) 



(16) 
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whereas the evolution of the orthogonal complement is not fully specified. The exception is 
for a spin-1/2 system. By unitarity, specifying one column of an 577(2) matrix necessarily 
constrains the other. Since our spinor lattice is described by a collection of noninteracting 
spin-1/2 subspaces labeled by quasimomentum g, if we specify a state-to-state mapping of a 
spinor Bloch state, we specify the SU (2) matrix on this block. We can achieve this using the 
state-mapping protocol defined in Sec. Ill that takes an initially localized Wannier state to a 
state distributed over a finite number of lattice sites. Such a map specifies a transformation 
on each Bloch state according to the Fourier relationship between the probability amplitudes 
in the Wannier and Bloch bases. Based on this relationship, we can use our state-to-state 
map to design a more general class of unitary maps on the wave function. 

To see this explicitly, consider the unitary evolution of an initial spin-up Wannier state 
(take I = without loss of generality), 

C/|0,|) = / dq(a(q)\q^)+P(q)\qA)) 

J -1/2 

oo 

= E (Q,tM)+Q4IU))- (17) 

l=— oo 

The quasimomentum functions a(q) and f3(q) in Eq. (16) are the Fourier sums of probability 
amplitudes in Wannier space, 

oo oo 

a(q) = e Q,t^ 27 r m = E w lMq - ( 18 ) 

l=—oc l=—oo 

As long as the Fourier transform of a(q) and j3(q) have support only over a finite extent in 
/, we can generate these functions by applying the state-mapping protocol of the previous 
section to synthesize the probability amplitudes q in Eq. (18). For unitary maps defined 
by a(q) and (3{q) whose Fourier expansion in Wannier states does not have a strictly finite 
support, more general control methods are required. 

We can easily generalize our result to include control through applied spatially uniform 
(possibly time-dependent) forces. In the TB approximation, expressed in the Wannier basis, 
the Hamiltonian for a linear gradient potential in dimensionless units takes the form 

oo 

H grad (t)= e F(t)L[i\n)(n\ + (i+si(t))\n){n\], m 

l=—oc 

where 5l(t) arises due to the off-set between spin-up and down lattices. We allow for mod- 
ulations of the overall force as studied in "shaken lattices" [38, 39, 40, 41, 42], and the 
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possibility of time-dependent variations in the relative positions of the two spin states, as 
could be implemented through modulations in the laser beams' polarization direction. The 
combination of this Hamiltonian, together with microwave-driven control described by Htb 
in Eq. (6), gives rise to a general unitary transformation that can be written using the 
interaction picture in the form, 

U(t) = D(t)U I (t), (20) 

where, 

/1/2 
dq\q-ri){q\®e- ixaz/2 (21) 
-1/2 

with X = Jo Sl(t')F{t')Ldt', rj = f* F{t')dt', and 



ijUi{t) = #/(*)£//(*), (22) 

where Hj{t) = D\t)HTB{t)D(t). The exact form of Hi is rather complicated, but all that 
really matters for our argument is that it is translationally invariant with the period of 
the lattice, TjHj(t)Tj = Hj(t). As a result, the solution, Ui(t) will be block diagonal in 
quasimomentum space, with the blocks consisting of SU{2) rotations. Thus by Eq. (20), 
the general unitary evolution will have the form, 

rl/2 

U(t) = e~^ 2 / dq\q-ri)(q\®U q . (23) 

J -1/2 

A control sequence of microwave driven rotations in a uniform lattice followed by a time 
dependent linear gradient can reach any unitary map of this form. 

We contrast this control with that achievable in a one dimensional sinusoidal optical 
lattice with time-dependent uniform forces and modulation of the lattice depth in the TB 
approximation, without symmetry breaking for right vs. left transport. Haroutyunyan and 
Nienhuis derived a general expression for the propagator in such a situation. [43]. In the 
quasimomentum basis, the map is 

U(t) \q) = e -^Wcos(2^-6(t)) | 5 _ ^ ( 24 ) 

where 

Jb(t) = f ji' od^M 1 ') 



a(t)e lb{t) = / dt'VL^y^) (25) 
Jo 

and Q(t) is the hopping rate between sites (symmetric to the left or to the right), and rj(t) 
is the same as above. The effect of the propagator is solely to induce a phase that varies 
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as the first order Fourier coefficient, in addition to shifting all of the quasimomentum by an 
amount rj. In contrast, Eq. (23) allows a broader class of unitaries to be synthesized. 

V. SUMMARY AND OUTLOOK 

We have described the control of transport of atoms through a microwave-dressed spinor 
optical lattice. Asymmetric lattices and spectral isolation provides a means to break the 
system up into a series of two-level systems, which aided the design of control routines. We 
restricted our attention to translationally invariant systems, with no local addressing, but 
the possibility of a uniform applied time-dependent force. Under these conditions, we can 
determine the constraints of reachable states and more general reachable unitary maps that 
take a localized atom at one site to a state coherently extended over n sites. Based on these 
constraints, we propose a constructive protocol for carrying out these control tasks through 
a sequence of SU(2) rotations acting in the two-level subspaces. 

An important consideration for practical implementation of our protocol is robustness of 
the control sequences to imperfections in the system. Because single particle transport is 
driven by a series of SU (2) maps, certain errors may be fixed by borrowing techniques from 
robust control of NMR systems[25, 26, 44]. In particular, spatial variations or miscalibrations 
in the microwave field strength and real or fictitious (i.e. light-shift induced) magnetic 
fields can be corrected with such techniques. On the other hand, the inevitable spatial 
inhomogeneities in the lattice potential can lead to spatial variations in the energy common 
to both of the levels, which causes an overall phase on the two-level systems that varies 
in space. Because this phase error is not an 577(2) map, it cannot be removed with the 
standard NMR composite pulse protocols nor their generalizations, and we must develop 
new methods to correct this in order to make our protocol robust to lattice inhomogeneity. 
This will be a topic of future investigation. 

There are a number of ways in which our method can be extended to more general con- 
trol of the transport of individual atoms. In the current work we restricted our attention 
to the single band, tight-binding approximation, in uniform lattices with two spin levels. 
More general protocols that do not restrict the bands can be used to study the control of 
coupled spin and spatial degrees of freedom in a broader context. Moreover, it has recently 
been shown [45] that the entire hyperfine manifold of magnetic sublevels is controllable with 
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applied rf and microwave waveforms, leaving open the possibility of combining the control of 
high dimensional spin and spatial degrees of freedom. In addition, breaking the translational 
invariance with quadratic or higher order potentials should allow the controllability of the 
system to be significantly enhanced. Other modifications of the spin-dependent potentials, 
such as the spatial variation in microwave transition frequency that arises in a strong mag- 
netic field gradient, can in principle allow spectral addressing of individual two-level systems 
and extend the controllability of the system. 

Finally, the techniques proposed here may also be extended to control the dynamics of 
many-body systems. For instance, it might be possible to use the microwave drive to synthe- 
size more arbitrary interactions between atoms than are dictated by the static Hamiltonian. 
Once interactions are included in the model, many-body unitary maps can by built from 
maps acting on restricted subspaces, as we have done here for single particles. Such tools can 
play an essential role for quantum simulations of many-body Hamiltonians, both for stud- 
ies of equilibrium properties such as the many-body phase diagram and non-equilibrium 
phenomena such as Lodschmidt echos [34] and the dynamics of phase transitions. 
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